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Abstract 

It has been shown that the cases of the JWKB formulae in 1-dim QM quantizing 
the energy levels exactly are results of essentially one global symmetry of both poten- 
tials and their corresponding Stokes graphs. Namely, this is the invariance of the latter 
on translations in the complex plain of the space variable i.e. the potentials and the 
Stokes graphs have to be periodic. A proliferation of turning points in the basic period 
strips (parallelograms) is another limitation for the exactness of the JWKB formulae. A 
systematic analyses of a single-well class of potentials satisfying suitable conditions has 
been performed. Only ten potentials (with one or two real parameters) quantized exactly 
by the JWKB formulae have been found all of them coinciding (or being equivalent to) 
with the well-known ones found previously. It was shown also that the exactness of the 
supersymmetric JWKB formulae is a consequence of the corresponding exactness of the 
conventional ones and vice versa. Because of the latter two exactly JWKB quantized 
potentials have been additionally established. These results show that the exact SUSY 
JWKB formulae choose the Comtet at al form of them independently of whether the 
supersymmetry is broken or not. A close relation between the shape invariance property 
of potentials considered and their meromorphic structure on the x-plane is also demon- 
strated. 
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1 Introduction 



It has long been known that for some number of potentials in 1-dim quantization prob- 
lems (and also in these cases of n-dim problems which can be reduced to 1-dim ones) their 
corresponding JWKB quantization formulae for energy levels (or some their generalizations 
Q or modifications ||) are exact |3|, ||, ||], whilst in most solvable cases (i.e. in those for 
which their corresponding energy spectrum is known by other means) the JWKB quantization 
appears to be only approximate. 

Evidently, the same solvable potentials which are quantized accurately by the JWKB 
formulae are also quantized exactly when the super symmetric (SUSY) modification of the 
JWKB method is used 0, [| ||]. But in general the SUSY JWKB quantization formulae 
similarly to the standard JWKB ones do not provide accurate quantization conditions for 
energy levels in most solvable cases of potentials [10|. 

Of course there have been attempts of proving that these known cases of both the con- 
ventional and SUSY JWKB quantization formulae |ll|] had to be accurate but from our 
point of view both these attempts suffered from arbitrary and erroneous conventions used 
to sum divergent series of necessary phases and one can easily convince oneself that some 
relevant parts of proofs in the thesiss mentioned are certainly incorrect (see also Sec.|3| and 
Appendix 1). 

A way to treat this problem properly is to use the Weierstrass infinite product represen- 
tations for mero- and holomorphic functions which the quantized potentials really are. But 
then it appears that it is necessary to take into account properly not only the phases coming 
from the infinite proliferations of (complex) turning points but also the phases coming from 
other (exponential) factors present in the corresponding Weierstrass products (see Appendix 

1)- 

Also a clear understanding of the exactness of the known (both conventional and super- 
symmetric) JWKB formulae seems to be still missing and they are considered to some extent 
as accidental. In particular, unknown are (possibly simple) criteria (different from a trivial 
direct comparison between JWKB results and exact energy spectrum known by other means) 
allowing us to conclude for which potentials the corresponding JWKB quantization formulae 
could be exact and for which certainly not. Having such criteria would be very important 
for applications of the JWKB method since it would relax us from appealing to other exact 
methods. 

The following are the aims of this thesis: 
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1. To prove rigorously the exactness of the known JWKB formulae; 

2. To clarify why the known JWKB quantization formulae are exact; 

3. To clarify the observed relation between the simultaneous (in-)accuracy of the JWKB 
and the SUSY JWKB quantization formulae; and 

4. To provide criteria allowing us to judge whether a given JWKB quantization formula 
can be accurate or not. 

For the beginning let us note however that in general there is no a common meaning of 
what a JWKB approximation really is. 

One of a typical way of considering it is just to get some standard semiclassical solution 
provided by the Schrodinger equation (SE) and valid in the considered domain. Next 
to truncate the corresponding infinite series on a given term. Such a truncation is then con- 
sidered as an approximate (to a given order in h) solution to SE in the domain chosen. In 
this way one gets (up to a ^-dependent constant) two types of the approximate semiclassical 
solutions (ASS) to SE correspondingly to two different signs of the classical momenta gener- 
ating them. Then correspondingly to the domain considered the dominate ASS is chosen. If 
none of the two ASS's is donaminating a linear combination of both is considered. 

Next, the ASS's constructed in the above way in different domains are smoothly joined and 
the ASS obtained in this way and covering the total domain of importance is then considered 
as a given order semiclassical solution to the problem considered. Taking the lowest order of 
this approximation we get what is then called the JWKB approximation to the problem. 

The above way of constructing of ASS has been described by Berry and Mount [^] 
for the 1-dim cases and by Maslov and Fedoriuk [^] for an arbitrary but finite number of 
dimensions. In particular the first paper discusses the difficult (and unsolved satisfactorily) 
problem of joining smoothly the ASS's across the so called Stokes lines on which they change 
their dominating character into a subdominant one (the so called 'connection problem'). 

From our point of view the above way of constructing ASS's suffers on a complete ig- 
norance of the exact solutions to SE which are approximated in this way and which in fact 
are unknown. This causes many confuses in applications of the method since together with 
the exact solutions also corresponding boundary conditions the solutions have to satisfy are 
ignored, too. 

Because of that and similarly to that it has been done in earlier papers [Jj], [D], 14, we 
shall be following here rather oppositely first starting with a particular set of well defined 
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solutions to SE satisfying a sufficient number of boundary conditions and having well denned 
semiclassical asymptotics (SA) as well. Since the latter property depends on domains the 
exact solutions are defined in we shall choose as such a set of them these having this property 
in a maximal way. 

It has been shown |l], [13], [L4|, [l5|] that such a set of the solutions to SE really exists and 
is known as a set of fundamental solutions (FS). In this thesis the descriptions 'ASS' and 
'JWKB approximation' shall be understood only just as the corresponding approximations 
to the fundamental solutions. This assumption has serious consequences for the form of the 
JWKB approximations which can differ seriously from the one described above. In particular, 
the presence of simple and second order poles in a considered potentials generate unavoidably 
changes in the corresponding JWKB formulae. 

A set of the FS's are accompanied by the so called Stokes graph (SG). Making use of 
the latter we get a uniform and systematic way of solving any interesting 1-dim problem 



both exactly and in the semiclassical limit |y, [L3|, 14, |15|]. The main property of SG is that 
it takes into account global features of a given problem considered in the complex planes of 
variables entered the problem (i.e. a position variable, energy, the Planck constant, some 
potential parameter(s), etc.). It is just these global features determining global structures of 
corresponding Stokes graphs which allows us to justify all the known cases of exact JWKB 
formulae as well as to get an insight what decides that a given JWKB formula can be exact 
or not. 

The fundamental solutions which play the main role in our approach have the following 
three basic properties: 



1. They satisfy some definite and desired boundary conditions [13, 14]; 



2. They possess well defined semiclassical limits [13, 14]; and 



3. The corresponding semiclassical series are Borel summed to the solutions themselves 
In fact FS's are the unique solutions to a given SE which have the above three properties 



altogether [15] 



The thesis is organized in the following way. 

In the next section we summarized essentialities related to Stokes graphs and fundamental 
solutions. 
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In Sec.^ we establish necessary symmetry conditions for potentials and their SG's to 
ensure the corresponding JWKB quantization formulae to be exact. Using these conditions 
we perform a systematic analyses showing that among all potentials which satisfy them only 
eight of them provide us with the exact JWKB quantization formulae. 

In Sec.|4| we describe the way of generalization of the results of the previous section invoking 
some our earlier results. 

In Sec.|| we argue that supersymmetric JWKB formulae cannot be exact if the conven- 
tional ones are not as such too and demonstrate why in all cases of the potentials of Sec.|| 
quantized exactly by the corresponding JWKB formulae the SUSY forms of the latter have 
to be also exact. By direct calculations we find in all these cases the validity of the SUSY 
JWKB quantization formula in the form of Comtet et al independent of that whether the 
considered superpotentials satisfy or break the supersymmetry conditions. The latter result 
is not however in a contradiction with that of Inomata et al [27] since our result concerns the 



exact JWKB quantization whilst that of the last authors is only the JWKB approximation. 



We discuss also in this section the result of Dutt et al [31| and Barclay et al [32] that 
the SUSY JWKB formulae are exact for the shape invariant potentials [ 30 1 and notice that 
all of them known as being shaped invariant under translational transformation are also 
quantized exactly by the JWKB formulae. According to that two more general theorems on 
the exactness of the SUSY and conventional JWKB quantizations of the shaped invariant 
potentials are formulated and proved. 

In SecJ^ we summarize our results and draw our conclusions. 



2 Stokes graphs, fundamental solutions and quantization 

We shall resume here the basic facts about Stokes graphs and fundamental solutions 



2.1 Stokes graphs 

Consider the SE written in the following form: 

E, A) - X 2 q{x, E, A)tf (s, E, A) = (2.1) 
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where: A 2mh 2 , q(x,E,X) = V{x, A) — E and a potential y(x,A) is assumed to be a 
meromorphic function of x and A with the following asymptotic behaviour for A — > +00 (h — ► 
0): 

y(x,A) ~y (x) + iyi(x) + ^y 2 (x) + ... (2.2) 



Together with q(x,E,A) we shall consider a function q(x,E,X) = q(x, E, A) + 5(x, E, X)/X 2 
where 5{x,E, A) behaves according to ( |2.2| ) when A — > +00. The precise form of 5(x,E,X) 
depends on types of singularities of q(x, E, A) in particular on whether the latter possesses 
simple or second order poles (see a discussion below). 

Let E be real and let Xi,X2, ■ ■ ■ , be roots of q(x, E, A) and yi,y2, ■ ■ ■ be its simple poles. 

Some of them can therefore be real but the rest ones are complex and conjugated pairwise. 

For each point Xi,yi,i = 1, 2, . . ., let us construct actions: 

X 

W[(x,E,X) = J ^/q(y,E,X)dy 

and (2.3) 

X 

Wf{x,E,X) = J ^/q(y,E,X)dy 

Vi 

and associate with them a system of lines defined by the conditions: 

MW[' p (x,E,X) = (2.4) 

These are Stokes lines (SL). To be a little bit more precise we call a Stokes line each 



connected set of points of the x-plane satisfying the conditions (^4). A collection of all 
Stokes lines is called a Stokes graph. 



If the ^-operation in ( |2.4D is substituted by the 9-one then the corresponding set of lines 
are called anti-Stokes lines (ASL). The two sets of lines are orthogonal to each other at all 
the points except the roots or poles of q(x, E, A). 
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Let Zk, k = 1, 2, . . ., be infinite points of the actions ( |2.3| ) i.e. the points where the integrals 
(|2.3j ) diverge to infinity. They are created by poles of q(x, E, A) (including these at infinities of 
the x-plane) and, therefore, their positions coincides with these of the poles. A total number 
of infinite points is assumed to be even infinite. 

The roots of q(x,E,X) and its simple poles are starting points for SL's. Due to our 
assumption about the root multiplicity only three Stokes lines can emanate from each Xi . On 
the other hand only one SL can emerge from each simple pole. Each of SL's starting from 
some root or simple pole: 1° can end at some other root or simple pole, 2° can end at the 
same root forming a loop around a second order pole, or 3° runs to an infinity point z^. 

A domain Sk containing a point Zk and bounded by some Stokes lines emenating from 
the roots of q(x, E, A) is called a sector of SG. Therefore, there are at least as many sectors 
as the infinite points zt- Typically the points zj~ collect a number of sectors which depends 
on a rate of increasing of the action when it approaches z^s. Sectors corresponding to finite 
Zfc's are also finite. The remaining ones extend to infinities of the x-plane. There are no roots 
Xi and no simple poles yi inside Sk, but there are some at its boundary. 

For the purposes of the thesis it is enough to know only the proper topology of relevant 
SG's i.e. their precise metric structures can be ignored. Therefore to draw the corresponding 
SG's it is sufficient for all the cases considered here to apply the following rules: 

1° From each root of q(x, E, A) (all roots are assumed to be simple) emanate three SL's and 
three ASL's; 

2° From each simple pole of q(x, E, A) emanates only one SL and only one ASL; 

3° From each second order pole of q(x, E, A) emanate only SL's or only ASL's depending on 
whether the pole coefficient is real negative or real positive, respectively; 

4° From each higher order pole (n > 2) of q(x, E, A) emanate n — 2 directions to which SL's 
are tangent asymptotically; 

5° Any SL and any ASL can have only single common point; 

The above rules together with the asymptotic properties of the actions ( |2.3j ) for x — > oo 
as well as an analytic behaviour of SG's on the parameters of q(x, E, A) the latter function 
can depend on allows us to draw any such SG considered in the thesis. 
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2.2 Fundamental solutions 



To any given SG we can attach to each of its sectors a solution to SE called a 
fundamental solution and having the following structure (2|, |l^, 14 , |l7| : 



* fc (x) = q-hx)e akXW * (x) Xk(x) (2.5) 



where: 



n>l 



2A 



2/1 2/n-l 

*/i J dy 2 ... J dy n uj{yi)uj{y2) ■ . . w(y n ) 

z k z k z k 



(l - e -*r k \(Wi(x)-Wi(vi))\ /j _ e -2«T fc A(Wi(»i)-Wi(«a))\ _ _ _ ^.6) 

^1 _ e -2<rkHW i (y n -i)-W i (y n ))'* 



with 



%) 1 g"(y) 5 «f 2 (y) 
= ~T7T ~ 7~ ITT + 7r~F77 ( 2 - 7 ) 



In the above formulae Xi is some of the roots lying at the boundary of Sk and o"fc = ±1 is 
chosen each time so as to ensure a sign of dt(ai~Wi(x)) to be negative for the whole sector Sk- 

One of the conditions determining the function 5(x,E,X) introduced earlier is to make 



all the multiple integrals in (2.6) convergent at their lower limits Zk- It appears that for the 
last reason this function has to be defined as non zero only if is a second order pole of 
the potential considered. The second reason appears when the FS is to be continued to a 
point being a simple or double pole for the potential. Namely, for both these cases we have 
to correct the potential always by the same term S(x, E, X) = (2(x — Zk))~ 2 at each simple or 
double pole of the potential. Of course, in the case of infinite number of these singularities 
the arising infinite series has to be sum into some function having them as its own simple and 
double poles. The 5-terms correcting the potentials considered in the above way we shall call 
the Langer corrections. 

We would like to stress at this moment that introducing the Langer corrections are un- 
avoidable part of the FS constructions whenever it is necessary to take into account the 
presence of the simple and double poles in the potential. 
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However, there are still another reasons for which particular forms of 5(x,E,X) have to 
be considered (see the next sections). 



If x £ Sk then all the integrations in ( |2.6| ) can be performed along so called canonical 
paths for which the condition $l(<JkWi(yj) — &kWi{yj + i)) < is fulfilled for any two successive 
integration variables. 



If x is any point such that the integrations in Q2.6| ) can be performed along some canonical 
paths then it is called a canonical point. A collection of all canonical points corresponding 
to the solution ^ is called a canonical domain. We denote the latter by D^. 

^fc's have the following two properties in their corresponding D^s: 
a. Their series (2.6) are uniformly convergent; 



b. Their asymptotic expansions when A — ► +00 are dominated by the first two factors in 

( |2.5| ) whilst the third one approaches then unity. 

Additionally we have: 

c. Every ^ is Borel summable in some Bk{Sk C Bk C Dfc) and ^fc's are the only solutions 

to SE with this property. 

The two first factors mentioned in the property b. above are the ones which just con- 
stitute our JWKB approximation to ^ we have talked about in the Introduction. But this 
approximation is valid only in the canonical domain of ^fc. 

All the fundamental solutions are pairwise independent. But since they are the solutions 
of the second order linear ODE ( |2.1| ) each three of them are linearly dependent. If for some of 
such a triad, say *$>j, canonical domains corresponding to them have common points 
pairwise then coefficients of a respective linear relation connected them can all be calculated 
in the following way @, |T§, : 

¥<(a?) = ai.^j(x)+a±^k(x) (2.8) 



where 



en, = lim - . , , ...etc 
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and x runs to Zk or Zj along canonical paths. The latter calculations allows us to get imme- 
diately corresponding JWKB approximations for the coefficients a^u^}.. 

2.3 Singularities of fundamental solutions 

Loci of singularities of the fundamental solutions in the x-plane coincide with the ones of 
q(x,E,X) and their nature is governed by the general rules (see for example [jf8|). Therefore 
if q(x, E, A) is holomorphic then such is each fundamental solution. The singularities at zeros 
of q(x, E, A) as provided by the representations ( |2.5| ) - (p7?| ) are therefore only apparent i.e 
they mutually cancel when the corresponding sums are performed. Nevertheless, they are 



real singularities in each factor in (2.5) as well as in each integral in (2.6) where they can 



cause troubles with taking the integrals (see the next section). 

Each simple pole y^ of q(x, E, A) is a source of a logarithmic branch point singularity for 
each fundamental solution which close to yk behaves as a(x — yk) + b(x — yk) \n{x — yk) with 
6/0. 



Each second order pole Zk of q(x,E,X) generates in each solution to SE ( |2.l[) a branch 
point of the form a(x — Zk) a + b{x — Zk) 13 if a — (3 is not an integer or the branch point of the 
form a(x — Zk) a + b(x — Zk)^ ln(x — Zk) in the opposite case. In the case of these fundamental 
solutions which have to vanish at Zk a real part of at least one of the numbers a, (3 has to 
be positive (one of the coefficients a and b can vanish in the case when the real part of the 
corresponding a or (3 is not positive). 

Each higher order pole Zk of q(x, E, A) generates a branch point which is simultaneously 
an essential singularity for each fundamental solution. The fundamental solutions which are 
defined in the sectors containing this Zk have to vanish at z^ in the corresponding sectors. 

It is easy to verify that the behaviour of the FS's (corrected if necessary by the Langer 
terms) around singular points of a potential considered is exactly such as described above 
sometimes being determined totally by their first two JWKB factors. However, let us stress 
it once again, that at the simple and second order poles this behaviour appears as a result of 
the Langer corrections. 

On all the figures of SG's drawn for the purposes of the thesis only cuts corresponding 
to the real branch points of FS's are marked whilst the branch points which follows from the 
F-F representation of FS's as given by (ET5|)-(|2.7|) are completely ignored. 



A set of all zeros Xk, simple poles yk and other poles z^ of q(x, E, A) determine for a real 
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E (and A) in a unique way a possibility of constructing a corresponding set of fundamental 
solutions among which there should be these two of them which determine the corresponding 
problem of quantization. If some (or both) of these two solutions cannot be constructed 
because of some specific properties of q(x, E, A) then the corresponding quantization problem 
cannot be formulated i.e. supposed bound states do not exist. The latter possibility is 
uniquely related to a particular pattern of SG which should be drawn in such cases. 



2.4 Quantization 



A quantization of 1-dim quantum systems with the help of the fundamental solutions has 
been described in many earlier papers flj, [H^, 14, 16 1. Here we sketch only the procedure for 
the case of two real turning points x\, X2 whilst the rest of them are complex and conjugated 
pairwise (we assume q(x, E, A) and E to be real). We assume also that our physical problem 
is limited to a segment Z\ < x < Z2 at the ends of which the potential has poles. In particular 
we can push any of Z\$. (or both of them) to =Foo respectively. 



To write the corresponding quantization condition for energy E and to handle simulta- 
neously the cases of second and higher order poles we assume Z\ to be the second order pole 
and Z2 to be the higher ones. 

It is also necessary to fix to some extent the closest environment of the real axis of the 
x-plane to draw a piece of SG sufficient to write the quantization condition. To this end 
we assume X3 and X3 as well as X4 and X4 to be another four turning points and 23 and Z3 
another two second order poles of V(x, A) closest the real axis. Then a possible piece of SG 
can look as in Fig.l. 
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Fig.l The SG corresponding to general quantization rule ( |2.10|) 



There is no a unique way of writing the quantization condition corresponding to the figure. 
Some possible three forms of this condition can be written as (l6j : 



exp 



A f q*(x,\,E))dx 



K 



Xi^4(A,£)x 2 ->3(A,£) 
Xl ^{\E)x2^(X,E) 



(2.10) 



and Xfc->j(A, E)k, j = 1,2,3,4 are the coefficients ( [E6| ) calculated for x — ► zj. The closed 
integration path K is shown in Fig.l. In the figure the paths 71— >3, 72^3) etc., are the 
integration paths in the formula (|2.6| ) whilst the wavy lines designate corresponding cuts of 
the :c-Riemann surface on which all the FS are defined. The same conventions in designations 
are maintained on the remaining figures 2-17. 



The above three forms ( 2.10 ) of the quantization condition are equivalent and can be 
substituted by another equivalent forms if the latter are still admitted by the SG considered. 
It means that in general we can choose between different forms of the conditions ( |2.10| ) 
according to our needs. In particular depending on the form of SG there are forms of ( 2.1 0| ) 
completely deprived of the JWKB phase factor on its LHS i.e. such forms are composed only 
from the x-f ac t° rs of FS's. 
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The condition ( p. 10 ) is exact. Its LHS has just the JWKB form. If we substitude each 



Xfc— >j(A, -E) in ( |2.10| ) by unity (which these coefficients approach when A — > +00) we obtain 
the well-known JWKB quantization rule. But in this way the latter is in general only an 
approximation to ( p.10 ). It is not in the following three cases only: 



1° All Xfc-»j(A, i£)'s in ( 2.1C| ) are really equal to (identical with) 1; 
2° They all cancel out mutually by some reasons; 

3° Both the above cases take place i.e. some of Xk-yj[\ E)'s satisfy 1° and some 2°. 



The first case is very rare and the only known example of it is just the harmonic oscillator 
potential (see below). It needs in fact for a given Xk^j(^, E) to have a possibility to deform 
its integration path properly to make all the integrations in ( |2.6| ) vanishing i.e. this condition 
demands some particular topology of turning points on the x-plane to happen. 

The next one if not happens accidentally can take place due to a possible reality of 
the coefficients entering the formula fl2.10| ) (where the coefficients can appear in pairs with 
their complex conjugations dividing them) or due to some possible symmetry of the potential 
V(x, A) relating the x _ coefficients present in the formulae. We shall show in the next sections 
that the latter case is the main reason for all the known cases of the JWKB formulae which 
provide us with the exact quantization conditions. In fact the symmetry properties of the 
potential as well as a particular topology of its turning point distribution cooperating together 
are the most frequent way of the realization of the JWKB formula exactness. 



3 Global symmetries of Stokes graphs and quantization 



The case when all Xj^k i n the condition ( 2.10| ) reduce to unity is exceptional and within 



the holomorphic potentials can contain only a single potential namely the harmonic one. 
(Note that for the holomorphic potentials we can always put 5(x, E, A) =0). 

To see this we note that for the case to happen it is necessary to have possibilities to 
deform the integration paths in the formula ( |2.6| ) for x's so as to make all the integrations 
vanishing. The latter property can happen if the integration paths can be pushed out to 
infinity i.e. none of roots of q(x,E,X) (which are branch points for the integrands in ( |2.6| )) 
can prevent such a deformation. It means that these roots cannot extend to infinity and that 
the corresponding SG for the case should look as in Fig. 2 with the blob containing all the 
roots of the case. Secondly, a number of sectors has to be limited to four as it is shown in Fig. 2 
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since only then the condition ( 2.1C| ) can contain only with vanishing integrations. Next, 
a number of roots inside the blob has to be finite since q(x, E, A) would vanish identically 
in the opposite case. Therefore, q(x,E,X) has to be polynomial. It is just a harmonic one 
since for the polynomial potentials a number of sectors exeeds by two a degree of a potential. 



Im x 




Fig. 2 The case of SG satisfied only by the harmonic 
oscillator potential among the holomorphic ones 



Therefore, other possibilities for the condition ( 2.10 ) to be the pure JWKB one are related 



to possibilities for x's in ( 2.10| ) to cancell mutually. Such cases to happen if not accidentall 



have to be related to some symmetries of corresponding SG's the latter being determind by 
relevant symmetries of underlying potentials. 

Suppose therefore q(x, E, A) to satisfy the following symmetry relation: 



q(y(x),E,\) = q(x,E,X) (3.1) 

for x — > y(x). We shall assume that it is always possible to find (if necessary) 5(x,E,X) 
such that ( |3.l[) is satisfied by q(x, E, A) as well under the same transformation. If the 
corresponding SG is also to be invariant under such a transformation then the full set of the 
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actions (|2.3|) which define this SG has to be invariant too, up to multiplicative constants. 



The latter freedom follows from the conditions (2.4) defining SL's. But according to ( |3.1| ) we 
have: 



x x y(x) 

J yJq(t,E,X)d£ = J Jq(y(0,\E)d£ = J yfq& E, A)^(£)d£ (3.2) 



from which we can conclude that the mentioned action set invariance is achieved if y'(x) = C, 
where C is real. 

Therefore the allowed transformations y(x) are linear. Since they constitute a group then 
it is easy to see that if \C\ 7^ 1 then q(x,E,X) has to have common accumulation points of 
their roots and poles. Because of that we shall limit our further considerations to less singular 
cases of q(x,E,X) what means that we shall put C = ±1. The latter limitation leaves us 
with only two types of the allowed symmetry transformations: the one which is essentially a 
reflection x — * —x and the other being a complex translation of the x -plane. 

Therefore the two resulting classes of potentials remaining invariant under the above two 
symmetry transformations are: a class of even potentials and a class of periodic potentials. 
Of course both the classes are not necessarily disjoint. It is, however, rather clear that the 



evenness of a potential alone is too week to ensure overall cancellations in (2.10) and it is 
just rather a periodicity of it which can work effectively to cause the relevant cancellations 
in (|2.10| ) to happen if it is possible at all. We shall show this below. 



3.1 Periodic holomorphic (entire) potentials 

In general q(x, E, A) as a meromorphic function of complex x can be periodic with at 



most two independent (in general complex) periods [19]. However, in the case of being 
holomorphic q(x, E, A) can have only one period (being a constant in the presence of the 
second one). Further, since q(x,E,X) is assumed to be real for its real arguments then its 
period can be only real or only pure imaginary. For the obvious reason we shall consider only 
the last case assuming for simplicity the period to be equal to 2iri. In this case q(x, E, A) can 



be expanded into the following Fourier series [19 



(x,E, A) = £ q n (E,X)e nx (3.3) 
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If the behaviour of q(x, E, A) at the x -infinity is to be of a finite type the series ( |3.3| ) has to 
be abbreviated providing us with a finite sum. The latter should contain at least three terms 
if we want q(x, E, A) to possess bound states. Let k and l(k > I + 1) be therefore the upper 
and lower limits of this abbreviation respectively. We shall consider just below in details a 
few cases of such abbreviated q's for which k — I = 2,3,4. By this we shall convince ourselves 
that the remaining cases of q(x, E, A) cannot provide us with examples of the exactly JWKB 
quantized potentials. 

In our investigations we shall make intensive use of the Weierstrass product representation 
for the abbreviated series ( |3.3] ) in order to perform necessary calculation of phases of q(x, E, A) 
alone as well as its functions. This represention is considered in Appendix 1. We have 
calculated there also explicitly, in order to provide us with an example of such calculations, 
the relevant total phases of q(x, E, A) for the case k = 2, 1 = considered just below. 

case: k = 2, 1 = 

We can write q(x, E, A) in this case as: 

q(x,E,X) = a(E,X)e 2x -2(3(E,\)e x + j(E,X) (3.4) 

where a{E, A), f3{E,\) and ~f(E, A) are known functions of E and A. In particular, for 
a = f3 = 1 and 7 = —E we get the well-known Morse potential p0| ]. 




Fig.3 The SG for the Morse type potential ( pip 



With a, j3, 7 > and f3 2 > 07 we get for q(x, E, A) = two real roots (modulo 2ni) and 
the corresponding SG shown in Fig.3 where x± = ln(/3 ± \J (/3 2 — 07) ). The quantization 
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condition ( 2. ICQ 



according to the figure looks now as follows: 



exp —\j)q^(x,X,E)dx 

K 



(3.5) 



It follows from the figure that X2^3 = X2^3 = 1 an d Xi^3 = Xi->3- The first of these 
identities is satisfied because both the paths 72^3 and X2^3 can De pushed out to infinities 
whilst the second because of the periodicity of the corresponding integrands in the formulae 
( |2.6| ) for xi^3 an< i Xi-»3- Therefore, we are left finally with the JWKB formula which gives 
exact energy levels in this case. 

It should be noticed, however, that the equality of the coefficients xi->3 and Xi->3 is n °t 
immediate i.e. it does not follow as a direct result of the periodicity of q(x,E,X). First we 
have to define the total phase of q(x, E, A) according to the prescriptions of Appendix 1 in 
order to define uniquely its square roots present in the coefficients X1-+3 an d Xi^3- For the 
case just considered it has been done in Appendix 1 where we have found that the phases of 
q(x, E, A) on the integration paths of Xi^3 an d Xi^S" differ exactly by 4-7T i.e. by the period 
of the square roots of q(x, E, A) just mentioned. 

It should be stressed at this moment that to get the last result the phases of the exponen- 
tial factor in the corresponding Weierstrass product (WP) have had to be taken into account 
i.e. counting the relevant phases providing by the roots of q(x, E, A) alone would give us 
incorrect result. 

case k = — I = 1 

In this case q(x, E, A) is given as: 



where a, (3, 7 all depend on E and A and are positive so that q(x, E, A) represents the infinite 
well. 



q(x, E, A) 



ae x + 7e 



—x 



2(3 



(3.6) 
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Fig. 4 The SG for the infinite potential well ( |3,G| ) 



For p 2 > «7(> 0) q(x, E, A) =0 has again two real roots (modulo 2ni) and the corre- 
sponding SG looks as in Fig. 4. For the corresponding quantization condition we can get the 
following two equivalent forms of it: 



exp 



A f qz (x, A, E)dx 



X2 ^ 3 (X,E) Xl ^- 4 (X,E) 



or 



(3.7) 



X2 ^ 3 (X,E) Xl ^ 4 (X,E) 



Xl ^ 3 (X,E) X2 ^ 4 (X,E) 



It follows from Fig.4 that Xl ^~ 3 (E, A) = X2 ^ 4 (E, A) = X2 ^(E, A) = 1, Xl ^ 4 (E, A) = Xl ^ 4 (E 
and X2 ^ 3 (E, A) = Xi-*i(E, A). The latter equality follows from the fact the phases of points 
of the corresponding integration paths differ by 2tt (by periodicity of q(x, E, A)) what causes 
the square roots of q(x, E, A) to differ by their signs on the paths compensated however by 
the opposite signatures of the coefficients considered. (Note that this time the corresponding 
phase difference calculations can take only into account the phase contributions of roots of 
(|3.6| ) alone since the exponential factor is now absent in WP corresponding to (|3,6|)), 

Taking all these into account as well as the following general equality p3| : 



X j->k(E,X) = Xk ^j(E,X) 
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we get finally for the quantization condition of the considered potential: 



exp 



-A J qz (x, A, E)dx 

K 



xUs(^E) (3.8) 



together with: |x2-+3(A,.E)| = 1. 



We have to conclude therefore that in this case the exact condition (3/7) cannot be reduced 
to the exact pure JWKB ones by the pure symmetry arguments only. 

The last potential, however, has been concluded initially by Rosenzweig and Krieger Q 
as being exactly JWKB quantized and next corrected by Krieger || as to be not. The main 
argument of the last author to support his conclusion was an observation of the unvanishing 
first order correction to the pure JWKB condition which follows from Q3.7| ). It is easy to see 
that in terms of the coefficient \2— >3{E, A) this argument means that the first integral in its 
representation (|2.6[) does not vanish i.e. the coefficient has to differ from unity. 



It is still worth to note that the potential (|3.q ) can be obtained from (|3,4| ) just by mul- 
tiplying the latter by e~ x , so that if the coefficients a, (3, 7 in both of them are chosen to be 
the same then their sets of roots coincide and their WP representations differ exactly by the 
above factor e~ x . And this is just this factor which introduces the dramatic difference be- 
tween the corresponding Stokes graphs of figures 3 and 4 and the corresponding quantization 
conditions. 

case k = 3, 1 = 
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Im x 




Fig.5 The SG for the potential ( |3.9| ) 

In order to satisfy the demand of two real turning points and the reality condition 
j(x, E, A) has to have now the following form : 



q(x,E,X) = a (e x -P + )(e x -p.)(e x + 1 ) 



(3.9) 



with a, /3±,7 > and all depending on E and A as usually. The corresponding SG is shown 
in Fig.5 where x± = ln(/3±) and x 7 = ln(— 7). The quantization condition is given therefore 
by: 



exp 



A f qz (x, E, X)dx 



K 



'xo^- 2 (E,X) X i^2(E,X) Xo-+z(E,\) 



(3.10) 



where the equality of Xi^2 and Xi-^2 to unity which follow from Fig.5 have already been 
used. Compairing further the phases of q(x, E, A) on the lines Qx = +7r and Qx = —it we 
find that they differ by Gtt. 



The condition (3.10) shows therefore that if there are not some accidental cancellations 
then the JWKB formula can be only an approximation in this case i.e. it cannot be exact. 



case k = 2, 1 = — 1 



This is the last case which we consider in details. The corresponding q(x, E, A) differs from 
pO|) by the last factor in which e x is to be substituted by e~ x . Thus q(x,E,X) represents 
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now an infinite potential well with the corresponding SG shown in Fig. 6. In comparison with 
Fig. 5 the horizontal SL emerging from x^ and its periodic distribution reverse their directions 
into the opposite ones. This change of SG does not however allow all the x' s to cancel in the 
corresponding quantization condition: 



exp 



A f (x, E, X)dx 



K 



(3.11) 



so that also in this case the JWKB formula cannot be exact. 

Im x 




Fig. 6 The SG for the infinite potential well: 

V(x, A) = a(e x - (3 + ){e x - 0-)(e~ x + 7 ) - a/3+/?_ 7 



corresponding to the quantization formula ( 3.11 ) 



Considering higher values of k — I (> 3) it is easy to note that all the conditions used 
as far (two real turning points, reality and periodicity) are not sufficient to course the full 
cancellations of the x' s m the corresponding quantization conditions. The main reason for 
that is that for higher values of k+l an increasing number of complex turning points in a basic 
strip of periodicity causes the x' s entering the corresponding quantization conditions not to 
be related any longer by the condition of periodicity, because all the relevant integrations are 
performed inside the same basic strip. 



The above situation does not change even if we make q(x, E, A) to be additionally an even 
function of x. 
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3.2 Aperiodic exactly JWKB-quantized potentials 

The periodic potential ( |3.4| ) which provides us with the exact JWKB quantization for- 
mulae ( |3.5| ) can serve also as the source of aperiodic exactly JWKB-quantized potentials. The 
latter can be obtained from the former by a trivial change-of-variable procedure x — * y(x) 
resulting in the following potential transformations: 

(3.12) 

x=x(y) 



q(x,E,X) 



q(x,E,X) 1 3 y" 2 (x) ly"'(x) 
y l2 (x) X 2 4y /4 (x) 2 y /3 (x) 



The only necessary demand for a relevant change is to provide by it in the resulting 
potentials a free (i.e. x-independent) term which can play a role of the energy parameter 
(see p8|], for example). 



The latter demand when applied to the potential (|3.4|) permits the following two possi- 
bilities: 1° e x l 2 — > x and 2°e x — > x. Adjusting properly a,/3 and 7 in (3.4) we get in this 



way the 1-dim harmonic oscillator potential and the radial parts of the 3-dim homogeneous 
harmonic oscillator potential in the first case and the Coulomb potential in the second one. 

The same method can be applied to the potentials which are not exactly JWKB quantized 
providing us with aperiodic potentials with the same property i.e. the method allows us not to 
make any further estimations of the resulting potentials for their being not quantized exactly 
by the corresponding JWKB formulae (what on their own could not be a simple task). 

Applying the method to the potential well ([3j]) for example only the substitution e x / 2 — ► x 
is allowed providing us with q(x,E,X) = ax~ 4 — (3x~ 2 — E with a,j3 > and — /3 2 /(4a) < 
E < 0. Of course, this potential can be considered as a radial part of a spherically symmetric 
potential with the repelling term x -4 and the attractive one x~ 2 the latter being affected by 
the centrifugal term contribution. It is seen that in this model a number of levels is finite 
being limited by the increasing of the angular momentum. According to our earlier result 
this potential cannot be JWKB exactly quantized. 

It is worthwhile to note that the above substitution when considered as the inverse trans- 
formation is nothing but the well known Langer change-of-variable procedure applied to the 
cases of the Coulomb and harmonic 3-dim potentials to obtain the exact JWKB formulae for 
the latter potentials. It is clear that in view of the above discussion Langer's substitutions 
are just what is necessary to be done in order to achieve the latter goal. 
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3.3 Periodic meromorphic potentials 

The reality condition demanded for q(x,E,X) allows us to choose for its two possible 
basic periods the one being pure real, and the second - pure imaginary. 

Within this class of potentials we can ignore obviously q(x, E, A) with a real period but 
without real poles. We have therefore to consider the following possibilities for q(x, E, A): 

a. It is holomorphic in some vicinity of the real axis but meromorphic outside it and being 

periodic with its unique imaginary period equal to 2iri; 

b. It is meromorphic on the real axis with the only imaginary period equal to 2m; 

c. It is meromorphic on the real axis with the only real period equal to 2ir; 

d. It is meromorphic on the real axis with two periods: a real one equal to 2tt and a pure 

imaginary one equal to iu with co being any positive real number. 



case a. 



In this case q(x, E, A) is assumed again to have (for some range of E) two real roots in 
its basic period strip defined by — ir < Qx < +tt. Let us assume also for a while that it has 
only four complex pairwise conjugated poles in this strip all lying on the imaginary axis (this 
position can always be achieved by a simple translation). If the order of these poles amounts 
to n then q(x, E, A) has to have the following form: 

q(x, E, A) = <&{x,E\) + q 2 (x, E, A) (3.13) 

< sinh 2 (x — ia) sinh ^ {% + io) \ 

where qi(x,E,X), i = 1,2, are holomorphic and periodic with their period equal to 2m and 
< a < it. Therefore roots of q(x, E, A) are given by the equation: 

f 1 1 } n 

gi(x, E, A) + j sinh -(x — ia) sinh -(x + ia) > -q2(x,E,X)=0 (3.14) 



with Qj's having forms of the abbreviated series ( |3.3j ). According to our earlier observations 
that more than two roots in the basic period streap prevent as a rule the JWKB formula to be 
exact we have to put in ( |3.14| ): n = 1, q±(x, E, X) = ae x + j3e~ x + 7 and q 2 (x, E, X) = const. 
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However, we have to consider the cases a ^ tt and a = it separately because their Stokes 
graphs differ in their structures. In particular, whilst in the first case q(x, E, A) will have 
simple poles at the points x = ±ia (mod 2-iri) in the second case the second order poles are 
generated. 

Satisfying finally the assumption of two real roots of q(x, E, A) we get the following two 
potentials: 



Vi{x) 



a\e x + Pi a\e x + P\ 



2 sinh \{x — ia) sinh \ (x + id) cosh x — cosh a 



V 2 {x) 



a 2 e x + P 2 
cosh 2 ^x 



(3.15) 



the second of which is essentially the Rosen- Morse one [plf| . 



To obtain from ( p. 15 ) the potentials which would have bound states some conditions on 



their parameters have to be satisfied. For the first of them they are: 



ai > > Pi 



(3.16) 



with the quantized energy E varying in the following range: 



2x(ai — x) 2 



< E < 



(|ai — x) — \y\) 2 + 2|y||ai — x\(l — cosa) 





(3.17) 



For the second potential in ( 3.15| ) 



we can put a 2 > > p 2 not loosing its generality so 



that the corresponding energy range is: 



01 



< E < 



(3.18) 



a 2 - p 2 
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Im x 




Fig. 7 The SG for the first of the potentials (|3.15|) 



A relevant Stokes graph corresponding to the first of the potentials ( 3.15| ) is shown in 
Fig. 7. It follows immediately from the figure that the JWKB formula: 



exp 



a±e x + ffi 
cosh x — cosh a 



E 



dx 



K 



(3.19) 



cannot be exact in this case since the coefficients Xi^4 an d X2^3 are n °t related by periodicity 
and do not cancel in the exact condition. 



The SG for the second of the potential ( 3.15| ) is shown on Fig. 8a. In order however to 
continue the relevant solutions corresponding to Sectors 1 and 2 to Sectors 3 and 3 (the latter 
two containing the second order poles at x = ±iri respectively) we have to choose properly 
the (5-piece of u> as defined by fl2.7|) to admit the integrals in (|2.6D to converge at the poles. 
One can easily convince oneself that the choice 5 = [4cosh(x/2)]~ 2 is sufficient for such a 
goal leaving simultaneously the original form of SG of Fig. 8a. unchanged i.e. it redefines 
the coefficient j3 into /?'(= — 1/(4A) 2 ) only. Then Fig. 8a provides us with the following 
quantization condition for the case: 



exp 



K 



\ 



a 2 e x + /3 2 



l 

uu 5 



cosh 2 tjX 



Edx 



' Xl-*z{KE) X 2^{\E) 



(3.20) 
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Fig.8 The SG's corresponding to the sec- 
ond of the potentials ( |3.15|) and the 
quantization formulae ( 3.20| ) (Fig. 8a) and 
O ) (Fig.Sb) 



According to Appendix 1 the total change of the phase of q(x, E, A) in fl3.20| ) is determined 
only by the distributions of zeros of its nominator as well as by the corresponding zeros of 
cosh 2 (x/2) being the denominator. Calculated (according to the rules of Appendix 1) with 
respect to the points of the lines = tt and Qx = —n (shifted by the period 2-7ri) the 
nominator phase change amounts to 2ir what is exactly the same as the total phase change 
of the denominator. Therefore, the total phase change of q(x, EX) is exactly equal to zero 
in this case what is enough for xi^3 and Xi^3 as wen as f° r X2^3 and X2^3 to coincide. It 
means that the RHS of ( |3.20| ) is equal to — 1 and the JWKB formula corresponding to ( |3.20| ) 
is exact in the case considered. 



It is interesting to note that the success of the JWKB formula corresponding to ( |3.20| ) to 
be exact depends completely on the total phase change of q(x, EX) to be equal to an integer 
multiple of 4tt. This condition permits also to produce the forms of the exact JWKB formulae 
different than the one discussed above. They can be obtained for example by doubling the 
number of turning points in the basic period strip. This can be done in many ways by 
choosing 5 properly. One of such choices is 5 = [4cosh(x/2)]~ 2 + a\2 sinh _2 (x/2) with real 
but arbitrary a / 0. This choice introduces to the basic period strip the second order pole 
at x = but also two additional zeros (lying close to the pole for small a). Both the zeros 
are real for a > or pure imaginary for a < 0. Let us note further that the pole at x = 
does not produce a singularity at this point for ^i(x) and ^2(x) but it can do it for their 
corresponding ^-factors (since it does it for the corresponding JWKB factors). However, 
the latter possibility does not affect the procedure of writing the corresponding quantization 
condition nor estimating periods corresponding to the coefficients Xi( x ) an d X2( x )- Choosing 
therefore for definiteness a to be positive we get the corresponding SG in the form shown in 
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Fig. 8b and the following quantization condition: 



exp 



K 



a 2 e x + (3 2 



16A 2 



cosh 2 |x 



sinh 



2 1, 



E 



dx 



X i^(\,E)x 2 ^(\,E) 



(3.21) 



Note that Xi-+3 an d X2^3 hi ( |3.21 ) are calculated on different sides of the (possible) cut 
emanating from the point x = 0. 



Now it is easy to see that the x's in the RHS of ( 3.2l| ) cancel mutually pairwise by 



periodicity and we obtain from ( p. 21 ) the exact JWKB quantization formula with a as an 



arbitrary real parameter. In particular we can put a = in the formula getting it again in 
the 'standard' Bailey's form Q. 

It follows from the above considerations that it is rather hopeless to look for other poten- 
tials of the case considered which could provide us with corresponding exact JWKB formulae. 
The proliferation of roots unavoidable for these potentials should prevent effectively the exact 
JWKB quantization conditions to appear closing the relevant quantizations inside the single 
period strip. 



case b. 



Assuming the presence of simple or second order poles in q{x, E, A) it is clear that we can 
allow only one such a pole in the main period strip. We assume its localization at x = 0. In 
this way the problem of quantization has to be reduced to a half of the real axis which we 
choose not loosing a generality to be the right one. The allowed classes of potentials should 
have therefore the forms: 



Vx(x,X) 
V 2 (x,X) 



a\e x + fa 
sinhx 

a 2 e x + Pi 
sinh 2 ^x 



(3.22) 
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If bound states are to exist for the first potential in (3.22) its parameters have to satisfy 
the following relations: 



01 < 0, -/3\<ai, ai + \l a\0l < E < 2a\ 



(3.23) 



V,(x) 



2a, 



a i+V a rP? 




"r V a rP 




Fig.9 The first of the potentials ( 3.22j ) and 
the corresponding SG 

A shape of the potential is shown in Fig. 9a. (Note that its local maximum is below its 
local minimum). The corresponding SG has to be modified because of the presence of simple 
poles in the potential and because of the corresponding quantization condition demanding 
for the wave function to vanish at x = 0. A 'minimal' choice for 5 is 5 = (2 sinh x) -2 so that 
the corresponding SG looks now as in Fig. 9b. 

The energy quantization demands now the solution *$>±(x) corresponding to Sector 1 and 
the solution ^fo(x)) from Sector to coincide. This provides us with the following quantization 
condition: 



exp 



A 



ai e x + Pt 1 



sinh x 4A 2 sinh 2 



E 



x 



dx 



Xl ^- 3 (E,\) X o->3(E,\) 
Xi_> 3 (£,A)xo-3CE,A) 



(3.24) 



From (3.24) it follows, however, that although the coefficients Xi-*3 an d Xi^s °f the 
formula mutually cancel (by periodicity) the remaining two coefficients do not (i.e. the latter 
coefficients are not real) and therefore the corresponding JWKB formula which follows of 
(3.24) cannot be exact. 
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a) 




V 2 (x) 










1 ) 


X 




X 




4a 2 










E 

?! 




















a 2+ p 2 











Im x 





1 Re x 



Fig. 10 The second of potentials (|3.22|) and 
the SG's corresponding to the quantiza- 



tion formulae ( |O0|) (Fig. 10a) and ( |3T27| ) 
(Fig.lOb) 



Not loosing a generality we can assume the parameters ct2 and /3 2 of the second potential 
in ( 3.22| ) to satisfy the following conditions: 



(3 2 , ol 2 + (3 2 > > 2a 2 + > ol 2 
< E < 4a 2 



01 



a-2 + 02 



(3.25) 



to ensure an existence of bound states in the local potential well. The shape of the po- 
tential is shown in Fig. 10a. The potential has to be modified by the 'standard' 5-term: 
5 = (4sinh(x/2))~ 2 to allow the construction of the FS at x = what results with the 
change: /3 2 — >• &i + (4A)~ 2 in the potential. The quantization condition corresponding to SG 
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of Fig. 10b reads now: 



exp 



K 



a 2 e x + j3 2 + 1 



TEX 1 



sinh 2 | 



E 



\dx 



(3.26) 



It follows from (3.26) and from Fig. 10b that in this case the coefficients Xi^2 and Xi^2 
cancel mutually by the periodicity arguments (the phase difference produced by the nominator 
of q(x, E, A) and equal to 2ir for the two integration paths 71^2 an d Ti^g i s cancelled by its 
denominator sinh 2 (x/2)) whilst the remaining two coefficients cancels by their reality (they 
are real and complex conjugated to each other). It means that the JWKB formula which 
follows from ( 3. 26] ) is exact. 



Again it is worth to note that the considered potential can be modified by the <5-function in 
a different way to generate at least four zeros in the basic period strip of the corresponding SG 
allowing the nominator of the modified q(x, EX) to change its phase by 4-7T between the earlier 
mentioned paths. This can achieved by putting 5 = (4sinh(x/2) -2 + aX 2 sinh _2 ((x — xo)/2)) 
with real and positive a but sufficiently small for xq satisfying x_ < xq < x + where x± are the 
outer two (of the total four) real turning points of the modified potential. The corresponding 
SG is then shown in Fig. 10c and for the quantization condition we get: 



exp 



A 



a 2 e x + (3 2 + 



wx 1 



sinh 2 § 



+ 



sinh 2 2^Ha 



E 



dx 



Xl ^- 2 (E,X) X 3^2(E,X) 
Xi^ 2 (E,X)x 3 ^ 2 (E,X) 



(3.27) 



It follows from ( |3.27D and from Fig. 10c that in this case the coefficients xi^2 and Xi^2 
again cancel mutually by the periodicity arguments whilst the remaining two again by their 
reality. It means that the JWKB formula which follows from ( 3.27| ) is again exact and 
coincides with the previous one when a = 0. 



The possibility of making the last modification enlarging the number of roots in the 
basic period strip to four still suggests to complete it differently, namely by adding the term 



coinciding exactly with the second of the potential ( 3.15 ). This of course needs also to add the 
corresponding standard 5-term to the potential obtained in this way. The potential we get in 
this way does not, however, satisfy the rule of no more than two turning points in the period 
strip so that a possibility of the exact JWKB quantization condition to appear should mostly 
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depend on symmetry properties of the relevant x~ coe ffici en ts. The considered potential can 
have bound states for the following regime of its parameters (see the formula below for the 
definition of the parameters): a, a' real and sufficiently close to zero and /3, (3' > 0. Then the 
SG corresponding to the case is shown in Fig. 11 and the quantization condition related to it 
is: 



exp 



A 



ae x + (3 + 



16A 7 



sinh 



2 x 



+ 



a'e x 



(3' 



i 

USA 2 



cosh 2 | 



E 



dx 



X1->2(E,X) X 3^2(E,\) 



(3.28) 




Fig. 11 The SG corresponding to the formula ( 3.29| ) 
quantizing the potential of Poschl and Teller 



It follows from the figure that the coefficients X3— >2 an d X3-+2 have to cancel mutually 
by periodicity but not the remaining two: none symmetry (except the complex conjugation) 
relates these two coefficients. However, when a = a' = the potential in (3.28) becomes 
invariant under the reflection x — > —x and then the coefficients Xi^2 and Xi^2 are equal just 
by the last symmetry. (Note, however, the role played in fulfilling this symmetry by 4-7T of 
difference between the arguments of q(x, E, A) corresponding to the case the latter takes on 
the paths 7i_>2 and 71^2)- 
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Therefore the following quantization condition: 



exp 



A 



+ 



1 

16A 2 



/?' + 



1 

16A 2 



sinh 2 | 



cosh 2 f 



rix 



(3.29) 



is exact. The potential in the above formula is of Poschl and Teller [22 



case c. 



The case contains the following four potentials: 



V 2 (x,X) 



ai sinx + (3 



cos x 

a 2 sin x + (3 2 
cos 2 £ 



. . . a'i sin x + /3{ a'/ cos x + 
y 3 (x, A) = ^ + -J : ^ 



COS X 



sinx 



— — j . . sin x + /3n a'o sin x + /3n 
F 4 (x,A) = — ^--^ + - 



cos 2 S 



sin 2 | 



7T 7T 

2 <S< 2 



-7T < X < 7T 



(3.30) 



7T 



< x < 



< X < 7T 



the second of which is essentially another of Poschl- Teller 




According to our earlier experience we cannot expect energy levels of the first potential 
(where ±«i + j3\ > 0) as well as of the third one to be exactly quantized with its correspond- 
ing JWKB formulae. This is because the points x = ±7r/2 which are singular for the first 
potential both lie inside its basic period strip and thererefore since the corresponding bound- 
ary conditions are formulated just for these points the ^"Coefficients entering the relevant 
quantization formula cannot mutually cancel i.e. the periodicity argument should not work 
in this as well as in the third cases. 

The second and the fourth potentials are more promising and as it is well known the case 
«2 = of the first one is quantized exactly by the JWKB formula. 

Consider therefore the first of them. In order to have the binding potential well we have 
to assume 02 > but the choice of sign of a 2 is arbitrary since both the cases are equivalent. 
So we shall put a 2 > for convenience. Next we have to notice however that asymmetry 
introduced to the potential by a 2 ^ completely eliminates the possibility of using the 
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periodicity arguments. Therefore we shall put «2 = in ( |3.30| ). Once more we have to 
choose 5 taking it in its 'standard' form 5 = (4cos(x/2)) -2 and getting the SG of Fig. 12a. It 
is seen from the figure that the coefficients x%—>-2 and Xi^2 as well as X1-+-2 anc ^ Xi—,2 are 
now equal by the periodicity arguments. Because of that the following JWKB quantization 
formula: 



exp 



cos 



2 x 



E 



dx 



(3.31) 



K 



is exact. 




Fig. 12 The SG's corresponding to the 
two variants ( ^3l|) (Fig. 12a) and ( ^32|) 
(Fig. 12b) of the JWKB formula for an- 



other Poshl-Teller potential 

Similarly to the hyperbolic cosine case there is again possibility to modify q(x, E, A) 
differently by putting 5 = a\ 2 sin _2 (x/2) + (4cos(x/2)) -2 . The first term of 5 (with an 
arbitrary real a) introduces two additional zeros in the basic period strip (necessary for the 
corresponding x's to have period 4-7r). The standard second term allows to construct the 
convergent solutions at the points x = ±7r. Choosing for definitness a > we obtain the 
effective form of the SG corresponding to these modifications as shown in Fig. 12b. It is seen 
from the figure that the coefficients Xi->-2 and Xi^2 as well as and xi^2 are equal 

by the periodicity arguments. Because of that the following JWKB quantization formula: 



exp 



cos 



2 x 



+ 



sin 



2 X 



E 



dx 



(3.32) 



K 



is exact for any a > and coincides with ( |3.31| ) for a = 0. 
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Considering finally the last of the potentials ( 3.30| ) first we remove (for the same reason 
discussed earlier) asymmetry in the latter putting a' 2 = cJ^ = and next we notice that in the 
basic period strip —ir < $tx < it the number of the four turning points is sufficient to make 
the relevant x' s periodic across the strip. Therefore the only necessary modification of the 
potential is the 'standard' choice for 5 i.e. 5 = (4sin(x/2)) -2 + (4cos(x/2))~ 2 what gives the 
SG shown in Fig. 13. The following relations come then out from the figure: Xi^3 = X— 1— >3 = 
Xi^3 and xo^3 = Xo^3- The first equalities in both of these equality sequences follows from 
the parity invariance of the potential considered whilst the second in the first one is satisfied 
by the periodicity arguments. Therefore the following JWKB quantization condition: 



exp 



K 



P' 2 + 



16A 3 



COS 



2 x 



+ 



T6A 2 " 



sin 



2 x 



E 



dx 



(3.33) 



is exact. 




Fig. 13 The SG corresponding to the exact JWKB formula ( |3,33j ) 



case d. 



Examples of the case are provided by elliptic functions The simplest candidates 

are the elliptic functions of the second order i.e. containing two simple poles or a pole of 
the second order in the corresponding basic parallelograms. In this case there are also two 
simple roots or a double one in each such a parallelogram of periods. However, in general, 
for a given elliptic function a number of its roots in each parallelogram of periods is always 
equal to the order of an elliptic function. Therefore, it is rather hopeless to expect an elliptic 
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functions with its order higher than two to be a good candidate of a potential to produce 
the corresponding exact JWKB quantization formulae. Considering, however, q(x, E, A) to 
be a second order elliptic function we have to assume that the two of its zeros are its real 
turning points x± placed between and 2ir (the latter being the real period of q(x,E,X)), 
i.e < x_ < x + < 2tt, whilst one of its poles is at x = 0. According to that we have the 
following two possibilities: 

a. there are two real simple poles of q(x,E,X): one at xo = and the second at x\, x + < 

x\ < 2tt; 

b. there is one double pole of q(x,E,X) (at x=0). 

Of course, in each of the above cases q(x, E, A) has to be completed to q(x, E, A) by the 
corresponding ^-function. The latter, however, has to be a sum of the second order elliptic 
functions with double poles at each singular point of q(x, E, A). Therefore S can be represented 
by a linear combination of two corresponding Weierstrass elliptic functions (case a) or should 
be proportional to such a function (case b). The periods of the Weierstrass functions coincide 
in each case with those of q(x, E, A). Therefore the SG's corresponding to the two cases have 
to look as in Fig.l4a,b respectively and the quantization conditions which can be prescribed 
to each of the graphs are the following: 




Xo^2(X,E) Xl ^2(X,E) 
Xo^2(A,£ki-2(A,£) 



(3.34) 



exp — X j> {a + -^^)p(x; 2-7T, iu) — E dx 

K 



X o^2(X,E) Xl ^- 3 (X,E) 
Xo ^- 3 (X,E)xi->2(\,E) 



(3.35) 
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Fig. 14 The SG's corresponding to the 
elliptic function F(x;2ir,iuj) (Fig. 14a) 
and Weierstrass one p(x; 2ir, iu) (Fig.l4b) 
taken as potentials 

Here F(x;2tt,iu;) is an elliptic function with two simple poles at xq and at x%, whilst 
p(x; 27r, iu) is the Weierstrass elliptic function. 



It follows from Fig. 14a that in the first of the above quantization formulae the periodicity 



arguments cannot work (x\ is not a shift of xq by 2n) nor in the form used in ( 3.34 ) nor in 
any other of its mutations. 



In the condition fl3.35|) one can use periodicity and reality arguments to show only that 
Xi^2Xo^3 = Xo^2Xi^3 ex P(— 2<5o->3) (where e>o^3 i s the phase of Xo^s) so that also in this 
case the RHS of this condition is not -1 i.e. the corresponding JWKB formula is not exact. 



4 More general exactly JWKB quantized potentials 

The periodic potentials considered in the previous section are the simpliest ones of the 
potentials quantized exactly by the JWKB formula. A generalization of their forms to the 
ones which still can keep the exactness of the corresponding JWKB formula can be done in 
the following way. 

Let V{x) means any exactly JWKB quantized periodic potential of the previous section. 
Let V(x,p) means a real parameter family of periodic (with respect to x) potentials with the 
property that in the limit p — > 0, V(x, p) approaches smoothly V{x) . Then, for p small enough 
the Stokes graph corresponding to V(x,p) has to resemble the Stokes graph corresponding 
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to V(x). By such a resembling we mean the following: 



1. To any singular point of SG of V{x) there correspond a set of singular points of SG of 
V(x,p) which reduce to the former point in the limit p — > 0. Each such a set we shall 
call a singular point blob (SPB); 

2. To any turning point of SG of V(x) there correspond a set of turning points of V(x,p) 
which reduce to the former point in the limit p —* 0. We shall call such a set a turning 
point blob (TPB); 

3. There is one to one correspondence between the sectors of the two Stokes graphs such 
that the sectors of SG of V(x,p) reduce smoothly to the corresponding sectors of SG of 
V(x) when p — > 0. In particular the boundary conditions are formulated for both the 
potentials V(x,p) and V(x) in sectors satisfying the correspondence just describe; 

4. Each set of Stokes lines which emerge from some SPB (TPB) can be mapped into a 
definite set of Stokes lines of SG corresponding to V{x) emerging from the point to 
which this SPB (TPB) reduces in the limit p — > 0. Each such a set can be divided 
into disjoint subsets of SL's each of them transforming smoothly when p —* into one 
particular SL emerging from the limiting point. 



Examples of V(x,p) with the properties 1 — 4 above with the Planck constant % as the 
parameter p can be found in Q. With these properties V(x,p) provides us with the JWKB 
formula quantizing exactly the energy levels of V(x,p). 



5 Some supersymmetric JWKB formula exactness 



In connection with the supersymmetric formulation of quantum mechanics the super- 
symmetric (SUSY) JWKB approximations have been suggested some of which being different 
from the conventional ones have appeared to be exact 0, H|. It has been also noticed however 
that their exactness have been parallel to the exactness of the conventional ones |7], ||, 11]. 



In the previous sections we have shown that the exactness of the conventional (i.e. not 
SUSY) JWKB formulae was rather exceptional and related to the periodicity properties of 
the corresponding SG's. Since the SUSY QM quantization problems seem to be governed by 
the same rules we can expect that the exactness of the SUSY JWKB formulae have to follow 
in some way from the traditional ones. We are going to show below that indeed this is the 
case. 



36 



Let us note, however, that there is also a common conviction that the SUSY JWKB exact 
quantization conditions are not only independent of the conventional ones but also their 
exactness in some cases of potentials is to be in contrast with the approximate character in 
these cases of the conventional JWKB formulae. As such potentials are considered the ones 
with the shape invariance property |3(| . We would like to argue below that also in these cases 
the parallelness of the exactness of both the kind of the formulae seems to be still maintained. 



Leaving the investigation of the latter relation for the later discussion let us examine first 
the question how the SUSY JWKB exact formulae follow from the conventional ones. 

For this goal let us remind that if a potential V(x, A) can be put in its SUSY form 
V(x, A) = V-(x\) = 4> 2 (x, A) — 4>'(x, A)/A + eo (eo is the energy of the fundamental level in 
V(x) if SUSY is exact) then the conventional JWKB quantization condition: 



-A § y / U(x,A) + %^ - Edx = (2m + l)m 
m = 0,1,2,... 



K 



(5.1) 



for the exact SUSY is to be substituted by jj], 



- A § ^/0 2 (x, A) — (E — eojdx = 2irim 

K 

m = 0,1,2, ... 



(5.2) 



If (|5.1| ) is exact then as we have mentioned (|5^) is very frequently also. 



Let us analyze how this can happen. We shall perform our analysis also for the cases of 
broken superpotentials <j) which can represent V(x, A). We shall show that in these cases the 
condition ( |5.2D remains unchanged if it is exact what is in contrast with its form representing 



the lowest JWKB approximation only in which case its RHS coincides rather with (5.1) 



At the beginning, let us note that because A can vary we can take it sufficiently large to 
expand the integrand in ( |5.2| ) into a series with respect to <f>' — 5/X. We get: 



Edx- 



n>l 



i r 


n 


i 

2 


\n-l 

A n!l 


r 

2 



K 



Ef 



—dx = (2m + l)m 



(5.3) 



where E = E 



eo- 
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Making further a change of variable: x — ► <j> = 4>(x, A) in the integrands of the series in 
flPD and putting A) = <f/(x(<t>, A), A) - 6{x((f>, A), A)/A and F 2 (</>, A) = <f/(x($, A), A) we 

obtain: 



-A 



Edx 



K 



n> 



i r 


n 


1" 

2 


\n-l 

A n!I 


l" 
2 



A', 4 



(02 -E) n ~h F 2 (<t>,\) 



(2m + l)vr (5.4) 



where the integrations under the sum in 



go now in the 0— plane. 



Next let us observe that for all the exactly JWKB quantized potentials considered above 
Fi2(<f>, A) are holomorphic functions of 4> outside some circles of a sufficiently large radius 
so that the circles contain the branch points at 4> = —y/~E and <f> = +y/E of the integrand 
denominators in (5.4). Moreover, both the functions grow with the same powers of 4> but 
not faster than the second ones. It follows then from ( |5.4j ) that all the integrands of the series 
in (|5.4j ) are also holomorphic outside such circles. The integrals can all be easily calculated 
then by taking the size of the contour K large enough and expanding their denominators: 



k>0 



-A 

K 

k + n — A 



kT 



n 



Edx-J2 



n>l 



i r 


n 


i 

2 


\n— 1 m 

A nil 


l" 
2 



A) d<j) 



02«+2fc-l F 2 (0,A) 



(2m + l)7n 



(5.5) 



A final result of the integrations in ( |5.5D depends now of course on the particular forms 
of the expansions of Fi^(4>, A) into their corresponding Laurent series. 



It can be easily checked that the series in LHS of (5^) becomes energy independent only 
in the case when the Laurent series expansions of F\i both abbreviate at least on the second 
power of 4>. This is just the case of the potentials considered. 



Suppose therefore that Fi j2 (0, A) = J2k>i ^i,2 ; fc(A)0 k + «i,2(A) + £>i,2(A)<£ + ci j2 (A)0 2 
Then from (^.5|) we get: 



K 



*P — Edx + iriSb o$c o + iri8 c o — 27ri — 

C2 



1-^-1 



(2m + l)vrz 



(5.6) 
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Below we shall do an inspection of the JWKB-quantization exact formulae of the previous 
section to calculate the corresponding coefficients a\ 2, &i,2 and c\ 2 as well as to show that all 
these formulae allow the quantization form (fT^) independently of whether the supersymmetry 
represented by <f> is exact or broken. This result, however, does not contradict the one 
obtained recently by Inomata et al (27| (see also pj| for further references) who have modified 
the Comtet et al formula ([|]^) with the aim to cover also the cases when supersymmetric 
potentials represent broken supersymmetry. They have argued that in such cases the 



RHS of (5.2) had to be transformed again into the conventional form of the RHS of (5.1). 
The source of the difference between both the conclusions is that our concerns the exact 
result whilst this of Inomata et al is only the lowest semiclassical approximation of the exact 
quantization condition. Nevertheless, as we shall see further that ([|]^) having the same form 
gives, however, different results for energy levels depending on whether the supersymmetry 
is exact or broken in the latter case reproducing effectively the result of Inomata et al . 

To follow further let us recapitulate all the potentials Vk (x) and the corresponding qk(x, E, X)- 
functions we have found in the previous section to be quantized exactly by the corresponding 
JWKB-formulae. They are: 



q~i(x,E,X) =Vx{x)-E = a 2 e 2x - 2(3e x - E, 



-00 < x < +00, (3 > > E 



q 2 (x,E,X) = V 2 (x) + 



4X 2 x 2 x x 2 

x,a, (3 > > E 



q 3 (x,E,X) =V 3 (x) + 



4X 2 x 2 



E, 



x- 



x,p,E > 



q±(x,E,X) = V±(x) 



(4A cosh I) 2 



E 
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ae 



i 

TEX 1 



cosh 2 § 



E, 



-oo < x < +00, P > 0, -j3 <2a 



q 5 (x,E,\) = V 5 (x) + 



(4Asinhf) 2 



E 



ae x + f3 + ^ 



TEX 1 



sinh 2 I 



-E, 



< x < +00, /3, a + /3>0>2a + /?>a, 



E, A) = F 6 (x) + (4Asi ^ hf)2 - (4A<x ^ f)2 ~E = 

ft + _ a + jix 7 _ E 

sinh 2 § cosh 2 § 

< x < +00, cc, (3 > 0, 



g 7 (z,£,A) = y 7 (x) + 



(4A cos I) 2 



— E 



a + 



TEX 7 



cos 



2 z 



-IT < X < tt, a > 



&(*, A) = y 8 (x) + IiX -L |F + IiX ^ r|F - S = 

or 1 9 «■ ^5 

cos^ 2 sir I 

0<x<7r, a, (3 > 



In order to represent the above potentials by their supersymmetric ones one has in prin- 
ciple to solve non uniform Riccati equations with their RHS given by the potentials listed. In 
general such a task is rather difficult. For most of the above potentials, however, it is possible 
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to find these representations just by a trivial guess. To each of the potentials listed above one 
can guess several (at least two) solutions one of which correspond to a superpotential realizing 
the supersymmetry exactly and the remaining ones corresponding to a broken supersymme- 
try. The latter means that the supersymmetry breaking can be realized in many ways. The 
ways considered below take into account only the possibility to define by a superpotential (f> 
the corresponding ground state solution by the following representation: 



^o(x) = exp 



A / 4>{y)dy 



a < x < b 



(5.7) 



where a,b (a < b) define boundaries of the corresponding quantization problem. Note that 
as given by (|5[ 



satisfies the SE (2.1) for E = €q with the potentials V(x, A)(= V-(x, A)) 
listed above. There are four possibilities: 



1 vanishes at both the boundaries a, b - the supersymmetry is exact and is the ground 
state wave function; 

2° and 3° \&o vanishes at one of the boundaries only (o or b respectively) - the supersymmetry 
has to be broken; and 

4° blows up at both the boundaries - the supersymmetry seems essentially to be broken 
but there is still possibility that the ground state has been constructed by the 
erroneous choice of — there are infinitely many solutions satisfying the SE considered 
with E = eo but blowing up at both the boundaries even if the corresponding ground 
state exists with this energy. 



The latter possibility cannot happen in the cases 2 and 3 : blowing up of a t one of 
the boundaries only means that the ground state with E = €q cannot exists in these cases. 
One can expect therefore that resulting relations between the energy spectra provided by the 
quantization conditions defined by the allowed superpotentials (f>k, corresponding to each of 
the potentials Vk, k = 1, . . . , 8, listed earlier, and the original spectra of the latter potentials 
can depend on the way the supersymmetry is broken by each particular <pk- 

Below we have enumerated all the allowed superpotentials 4>k corresponding to each of 
the potentials Vk, k = 1, . . . , 8, with the properties 1° — 4° just discussed (attaching to each 
of them the corresponding category) together with their iq^-functions: 
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4>i(x, A) = \a\ 

Wi) = Wi) 

h = 1; e = 




A) = -(ale 21 + 



Al- 
ia] 2A' 



Fi(0i) = F 2 (0 1 ) = 1 -^ + ^-, 

|a| zA 

61 = 15 eo = - ( rr + 7x) 2; 

a 2A 



, , |2Z + 1| + 1 

F 1 (0 2 ) = A(2|2/ + 1| + 1) 



Aa 

|2Z + 1| + 1' 

aA \2 
2 - |2Z+l|+lJ 



(|2Z + 1| + 1) 2 ' 

aA \2 



F 2 (0 2 ) = 2A- 



2/- 



-l\+l> 



. 2|2Z + 1| + 1 
c l = A 7TTTT— TT— TT2 1 c 2 



|2Z + 1| + 1 ' 

2A 



(|2/ + 1| + 1) 2 ' ^ |2Z + 1| + 1' 
(Aq) 2 



eo 



(|2Z + 1| + 1) 2! 



1(1 + 1) 
A 2 ' 



Z < -1, I > 0; 



</> 2 (x, A) 



|2Z_+lj -1 
2Ax 



Aa 

|2Z + 1| -1' 



Fi(^) = -A(2|2Z + l|-l)- 



i>2 + 



aA 
|2Z+1|- 



if 



F 2 (<h) 



-2A- 



t>2 + 



(|2Z + 1|-1) 2 ' 



aA 



|2J+1|-1- 



|2Z + 1|-1 ' 
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Cl 



. 2|2Z + 1| -1 



2A 



(|2i + l|-l)2' |2Z + 1|-1' 

(Aa) 2 



1(1 + 1) 
A 2 ' 



(|2Z + 1|-1) 2 ' 

I <-l, Z > 0; 



1° 



Fi(0 3 ) = A(2|2Z + 1| + 1) 
F 2 (<fa) = A 



4>3(x, A) = \a\x — 
|2Z + 1| + 1 



|2Z + 1| + 1 



2Ax 



+ 2\a 



A 



2 03 



^ 2 + 2H^i)i 



+ 



a 



6 3 + 2 M ^ 



2(|2Z + 1| + 1) 2 ' 2|2Z + l| + 2' 

§fr + (^ + 2|a|J*±^)* 



o 2|2Z + 1| + 1 o 
Cl - A (|2Z + l| + l) 2 ' ° 2 



|2Z + 1| + 1 

2A 



eo = (|2Z + l|+2)M, f3= l -V±V 



|2Z + 1| + 1' 

Z < -1, Z>0; 



2° we get the case from 1° substituting there |2Z + 1| by — |2Z + 1|; 
3° we get the case from 1° by the substitution \a\ — > — \a\; 

4° we get the case from 1° substituting there |2Z + 1| by — |2Z + 1| and \a\ by —\a 



1° 



|2/ + 1| - 1 x 
h = 1 7T tanh - + 



4Aa 



4A 



F 1 (</» 4 ) = -A(2|2Z + l|-l) 



2 |2Z + 1| -1' 

i 

4\a \2 



I2Z+1I-1- 



|2Z + 1|-1 



+ 



2|2/ + 1| - 1 



16A 



FM ~ " 2A |2Z + 1|-1 + 8A : 

2|2Z + 1| -1 2A 
C1 -" A (|2Z + 1|-1) 2 ' C2 ~ 

"|2Z + 1|-1 2Xa 



|2Z + 1| -1' 



eo = - 



a + (3 



2A 

1(1 + 1) 
(2A) 2 ' 



|2Z + 1| -1. 

Z < -1, Z;Z > 0; 
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4° we get the case from 1° by the substitution |2Z + 1| — ► — |2Z + 1| 



1° 



|2Z + 1| + 1 iU x 
■ coth — 



4Aa 



Fi(0 5 ) = A(2|2i + l| + l) 
F 2 (0 5 ) = +2A 



4A 2 |2Z + 1| + 1' 

05 + p^ppr) 2 2|2Z + 1| + 1 



(|2Z + 1| + 1) 2 16A 
(05 + p^r) 2 |2Z + 1| + 1 



|2/ + 1| + 1 
. 2|2Z + 1| + 1 



8A 



2A 



(|2Z + 1| + 1) 2 ' z |2Z + 1| + 1' 



a|; e = - 

1(1 + 1) 



|2Z + 1| + 1 2Aa 



+ 



a + /3 



(2A)2 ' 



2A ' |2Z + 1| + 1 
1<-1, l>0; 



4° we get the case from 1° by the substitution |2Z + 1| -> -|2Z + 1| 



1° 



[2/ + 1J-1 x |2T + 1| + 1 x 
66 = 4A tanh 2 4A C ° th 2 ' 



Fi(^) = --(2|2Z+l|-l) 



^ + (2|2/ + 1|-1) 



Igr+jjj-i 

(2A)2 



5 06 + Ul + (2|2i + 11 — 1) 
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6| + (2|2I + 1|-1) 
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(2A) 2 



0l + (|2Z + l|-l)Wi 



^2(06) = -A 
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c 2 — 



(|2Z' + 1| + 1) 2 

2|2Z + 1| - 1 
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2A 



|2Z + 1| -1' 



-+ 
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_(/-/'- l) 2 _ 1(1 + 1) _ /'(/' + 1) 
60 (2A) 2 ' ° (2A) 2 ' ^ (2A) 2 ' 

|2Z + 1| - |2/' + l| > 2, M'<-1, /,/'>0; 

2° we get the case from 1° taking Z, /' satisfying \2l + 1| — |2/' + 1| < 2 or substituting |2Z + 1| by 
-\2l + 1| there; 

3° we get the case from 1° substituting there \2l' + 1| by — \2l' + 1| and next taking /, /' satisfying 
|2i' + l|±|2Z + l| >2; 

4° we get the case from 1° substituting there \2l' + 1| by — \2l' + 1| and next taking /, /' satisfying 
\2V + 1| ± \2l + 1| < 2; 



1° 



12/ + 11-1 x 
h = Tx tan-, 

F l( *) = -A(2|21 + l|-lV # 2|2/ + 11 " 1 



F 2 (<h) = -2A 



(|2/ + 1|-1) 2 16A 
6? \2l + 1| - 1 



|2Z + 1|-1 8A 

2|2/ + 1|-1 2A 
c l = ~^ 7T7T1 , 71 7777 > c 2 = ~ 



(|2/ + 1|-1) 2 ' z |2Z + 1|-1 
|2/ + 1| -l x2 



eo 



4A 



4° we get the case from 1° substituting there |2/ + 1| by — 12/ + 1|; 



l u 

12/ + 11 + 1 x 12/' + 11-1 x 
— tan — cot -. 

4A 2 4A 2' 



= -(2|2/+1|+1) 



2 + (|2/' + l|-l) 



|2/ + 1| + 1 
(4A) 2 



+ 



i> 2 + (i2/'+i|-i)^j+i; 



~(2|2Z' + l|-3) 



,2,^ + 1) 

3 s ~r 



|0 8 -U 2 + (|2/' + l|-l)^±i 



(|2/ + 1| + 1) 2 
*+i|+ 



"2{<PS 



= A 



0l + (|2/' + l|-l) 



A 2 

|2/ + 1| + 1 

(4A) 2 



+ 



(|2/' + l|-l) 2 

0l + (| 2 r + i|-i) 



|2i+l|+l 



|2/ + 1| + 1 
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0| + (2|z' + i|-r 



12/ + 11 + 1 



(4A)= 



6§ + (|2Z' + l|-l)l^ 



A 



41' - 1 



(21 



A2 



4/ + 3 
'(2/ + 2) 2 ' 



f 2 



(|2Z' + 1 
A 



1) 



c 2 



A 



eo 



(l-l'+lf 
(2A)2 



Q: 



Z(Z + 1) 
(2A)2 



P 



V{V + 1) 
(2A) 2 



Z + l' 

M'<-1, Z,Z' >0; 



2° we get the case from 1° substituting there |2Z + 1| by — \2l + 1|; 
3° we get the case from 1° substituting there \2l' + 1| by — \2l' + 1|; 

4° we get the case from 1° substituting there \2l + 1| by — \2l + 1| as well as \2l' + 1| by — \2l' + 1|. 

In the above calculations only the non vanishing coefficients a% t 2, 61,2 an d c+2 at the 
highest power of <j) have been given. Now we can use them to calculate the three pieces of 



the LHS of (5.6) and to convince ourselves that in all the cases considered all the pieces 



contribute the total value m\ only so that (|5.5|) always reduces to (5.2) 



A comment to these calculations is necessary in the cases of the V3, Vq and Vs potentials. 
The necessary integrations that lead to ( |5.6j ) are performed here on the two sheeted Rie- 
mann surfaces of the variable on which the corresponding functions F\^{4>) have different 
asymptotic properties for — > 00. The superscripts '0' and '00' at the coefficients ai.2, 61,2 
and ci 2 indicate that the integrations contributed to them have been performed on the two 



different sheets. In such cases the coefficients with both supercripts contribute to (5.6) but 
these contribution should be multiplied by 1/2 each (see Appendix 2 for an example of such 
calculations) . 



Let us finally note that compairing the energy levels obtained by the formula (|5.l| ) with 
those obtained by ( |5.2j ) using in the latter the respective superpotentials of the cases l°-4° 
we get the result that the energy levels given by (5.1) are reproduced by (5.2) 



i. exactly in the case 1° of the superpotentials; 

ii. being shifted up by a half of a unit used to enumerate the levels in the cases 2° — 3° of the 

superpotentials; 

iii. being shifted up by a whole unit used to enumerate the levels in the cases 4° of the 

superpotentials; 
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It is clear that the above differences follow as a result of the different enumeration of 
energy levels in the compaired spectra (m in (fO]) starts from zero whilst in ( |5.2|) from unity) 
as well as due to different choices of the energy levels eo with respect to which the levels of 
the spectra are measured in every of the cases l°-4°. 



5.1 SUSY and conventional JWKB quantization of shape invariant poten- 
tials 



We have shown in the previous section that all the exactly JWKB quantized cases of 
potentials are also quantized exactly by their SJWKB quantization formulae. However, the 
latter property of the considered potentials has been established also as being closely related to 
their common property of being shape invariant |3(J]. The latter means that each Vk(x, A) = 
VL ( fc(a?, A), k = 1,...,8, depends additionally on some parameter a so that for its SUSY 
partner V + ^{x, A, a) we have: 



14 



V- k (x, A, a) + R k (ai) 



(5- 



k = 1, 



with ai = /fc(a). In the case of the considered potentials each is simply a translation of 
the parameter a. 

The exactness of (|5.2|) following from (pTsj) has been suggested by Dutt et al pi] and 



established by Barclay et al [Q. It was argued also (see Cooper et al and [29], for example) 
that the exactness of SJWKB formulae (|5.2| ) takes place even when the exactness of the 
conventional ones fails. The latter claim, however, needs not be necesserilly true and we 
would like to show below that in the case of the translational shape invariance all the known 
cases of the potentials are JWKB quantized exactly, too. To this aim let us note that on the 
list of the eight of them cited above there are two still lacking on the list when compaired 
with the corresponding list of Cooper et al pS|], These two are: 



. , . a + 3 sin x ir n 

V 9 (x, A) = = , -- < x < +-, a > 3 > 

cos^ x 2 2 



V 9 mm =zW^-8i + a) (5.9) 
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a + /?sinhx 

Vio(£, A) = g > — 00 < x < +°°) p > 

cosh x 



One can easily convince oneself, however, that completed by the 'standard' <5-terms 
((2 cos a;) -2 for the first potential and —(2 cosh x)~ 2 for the second one) the potentials are 
exactly JWKB quantized no doubt reflecting the fact that this exactness follows in some 
although not obvious way from the underlying periodicity of their SG's shown in Fig. 15. 




Rex=-x 2 Rex=x. Rex=-x, Rex=x 2 

Fig. 15 The SG's corresponding to the po- 
tentials Vg(x, A) (Fig.l5a) and Vxo(x, A) 
(Fig. 15b) given by the formulae (]5[ 



e x=-jt-x, He x=-jt-x, 



Imx 



o , O 



.:::<:=>:. Kcx=rt-x, Hex=Tr-x, 



Rex 



This exactness can be, however, concluded also applying back the procedure described in 
the previous section i.e. the exact JWKB formulae for the considered potentials follow from 
the exactness of the corresponding SJWKB ones. To see this let us find the superpotentials 
09,10 and their Fi^-functions corresponding to the potentials ( |5.9[ ), They are: 



|2Z + l| + |2Z' + l|+2 12/ + II - 12/' + 11 

t>9 = — tan x + 



4A 



4A cos x 



bl + a 2 ) 



+ a 2 )a' - ah 1 + bfaJ<t>§ + 



b 2 



F 2 (<j> 9 ) = (4 + a 2 ) 



(b 2 + a 2 )(f>l + a 2 (a 2 - b 2 ) + 2a&0 9 yVl + a 2 - b 2 
((pi + a 2 )a - ab 2 + b(j) 9 ^l 



b 2 + a 2 - b 2 



„°°i 



(b 2 + a 2 )cp 2 + a 2 (a 2 - b 2 ) + 2ab(j) 9 ^ 2 + a 2 - b 2 
111 _ 1 1 1 



„0O2 



a + b 4A (a + b) 



2 ' °i 



a - b 4A (a - b) 2 ' 
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~°oi __ \ 002 



a + b' 2 a-b 

« + P = ^ ' «-P = ^2 ' 

12/ + II + \2l' + 1 1 . 1 

+ 7TT' a = a - — , 



4A 2A' 4A' 

|2Z + 1| -\2l' + l\ , . , . 

— , u < -i, m >o 



4A 

2° we get the case substituting in 1° \2l + 1| by — \2l + 1|; 
3° we get the case substituting in 1° \2l' + 1| by — \2l' + 1|; 

4° we get the case substituting in 1° \2l + 1| by -\2l + 1| and \2l' + 1| by -\2l' + 1|; 



I2Z + 1I-1 , b 
>io = tanh x + 



2A cosh x 



- a 2 ) 2 + 4A(0? O - a 2 ) fa(^ - a 2 - 6 2 ) - ib^ <t> 2 w - a 2 - b 2 

Fl( ^ = it ; / % 

4A (ft 2 - a 2 ){(/) 2 - a 2 ) + 2a 2 b 2 + 2iab<t> W yJ<f> 2 10 - a 2 - b 2 



F 2 (<p 10 ) = (4>i -a 2 ) 



, 2 a(0io - a 2 - b 2 ) - ib(f> W \l 4> 2 W -a 2 - b 2 



{b 2 - a 2 )(^ 2 - a 2 ) + 2a 2 6 2 + 2iab<f> 10 J<f> 2 10 - a 2 - b 2 



a-*6 4A(a-*6) 2 ' 1 a + i& 4A(a + i6) 2! 



(a — ib' a + ib 

b>0, l>0, 1<-1; 

4° we get the case substituting in 1° \2l + 1| by — \2l + 1| or allowing I to vary in 1° in the segment 
— 1 < I < (the allowed b is then negative in both the cases). 

Since the behaviour of the functions in both of the above cases corresponds exactly to 
the assumptions we have done about them earlier to relate the conventional and supersymetric 
JWKB formulae these relations have to be maintained also in both the cases considered. A 
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detailed demonstration of this is performed in Appendix 2. The latter illustrates a typical 
way of getting such relations. 

All the results we have obtained in this section are the good illustrations of theorems we 
are going to formulate and to prove below. 

Let V(x, A) = V-(x, A) satisfies the following assumptions: 

1° V—(x, A) being shape invariant is meromorphic on the x-plane having there finite second 
order poles and diverging at infinity quadratically or exponentially with x (the two 
latter possibilities excludes each other of course); 

2° The map: x —* <p = 4>{x) defines a finitely sheeted Riemann surface of the c/>-variable; 

3° All turning points of q(x) = V±(x,X) + — E are transformed into pairwise i.e. 
one pair into one sheet of R^ and each pair can be used to define equivalently the 
contour K in the quantization condition ( |5.1| ) (that is each such a contour surrounds 
the corresponding turning point pair); 

4° The functions -Fi(</>) and F%{<j)) defined on R^ are holomorphic on R^ outside some circle 
of sufficiently large radius not branching at infinities of any sheet. 

The above assumptions allow us for the following conclusions: 

a. 4>{x) maps the contours K of the assumption 3° into respective contours on R^ which 

surround the corresponding maps of the turning point pairs on R^; 

b. R,)) is cut with the branch points 4>{ x k) satisfying: <p'(xk) = where x/- are finite regular 

points of 4>(x). All these points lie outside all the contours Ka, and are the square root 
branch points for F\ and i 7 ^- (The latter type of branching follows from the equalities: 

if 4>'(xk) = then close to x^ (f>'(x) ~ a(x — Xk) and 4> ~ ol{x — Xfc) 2 /2); 

c. If (j>'4 vanishes on R^ at some of its regular point 4>q linearly i.e. (f)'(x) ~ a((f> — 4>q) then 

the point is a map of an essential singularity of (j>(x) lying at infinity i.e. <p approaches 
4>o exponentially: 4>{x) m <J)q + Ce a ' x ; 

d. Second order poles of V±(x, A) as well as their infinite singular ones are transformed into 

infinities of different sheets of R^. The divergence to infinity of the F\ ^-functions is 
the following: 

i. For V-(x,X) diverging as e x for x — > oo Fi% diverge linearly with <f) when <fi — > oo 
on a given sheet; 
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ii. For V-(x, A) diverging as x 2 for x —> oo iq^ approaches constant values when 



4> — > oo on a sheet which vicinity of its infinity is a map of the corresponding 
vicinity of the pole xq; 

Let us add yet to the four above the following one more assumption: 
5° The integrand of the following integrals: 



where K$ is any of the contours of the remark a. above do not possess outside the contours 
singularities different than those described in the conclusion b.. 

From the assumptions 1° — 5° and from the remarks a.-d. the following two theorems 
come out: 

Theorem 1 The SJWKB formulae with the superpotentials <f>(x, A) corresponding to V±(x, A) 
are exact independently of whether the supersymmetry is exact or broken. 

Theorem 2 The conventional JWKB formulae for V±(x, A) are exact. 
Proof of the Theorem 1 . 

The theorem follows from the following sequence of equalities: 



iii. For x close to a second order pole of V-[x, A) iq ; 2 diverge to infinity as <fi- 



> 2 when 




(5.10) 




K 



K 



+ <j>(f(Fi(ai), E - R( ai )) - f(F+(a),E))dx = ... 



K 




(5.11) 



K 
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+ jr<f(f(F 1 -(a p ),E-R(a 1 ) 

P =1 K 



R(a p ))- 



/(F+(a p _i), E - R(ai) R(a p ^)))dx, 



where f(F 1 ,E) is defined by: 



ao = a, R(ao) = 



A' 



4>\a) ± l^'(a) + 



A 2 



(a)-E)*dx+ f f(F?(a),E)dx (5.12) 



if 



A' 



with = ±0' + «J/A. 

From assumption 3° it follows that every of the contour integrals in the sum of the RHS 
of ( |5.11| ) when rewritten to be taken on some sheet of R^ can be taken on each sheet of R^ 
in the following way: 



f(F ± (x,a),E)dx 



A 



1 n 

-E 



A,,, 



(5.13) 



where ^(0) = ^' {<ft)) ■ 



Now it follows further from assumption 5° that every contour K^ r 



, n, can 



be deformed on a sheet which it is defined on to a circle of sufficiently large radius and to 
pieces of this contour which cancel mutually with analogous pieces of other contours. The 
net result of these deformations are the integrations performed on every sheet along the 
circle with sufficiently large radius. Outside the circle the integrated /'s (divided by F2) are 
holomorphic and diverging to infinity not faster than the second power of (p. This guarantees 
that all these integrals can be calculated in the way similar to that used by us earlier. It is 
easy to check that independently of the type of the divergencies listed in the points i.-iii. 
above each infinity contributes the same to the sum ( 5.13] ) namely ^in/X for the F^-cases 
respectively. Therefore, the total value of the integral in the LHS of (|5.13|) is also ^iir/X 
accordingly. Finally the formula ( |5.11 ) becomes: 



: (a) — E)2dx 



E + R(ai) + . . . + R(a m ))2dx + 2mm 



A 



A 



ao = a, R(ao) = 



(5.14) 
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Putting now in ( |5.14|) E = R(a\) + . . . + R(a m ) = E m we get the result Q5.2| ) where for 
the broken supersymmetry the integer m starts rather from m = 1. QED. 



Proof of Theorem 2. 



The claim of the theorem follows immediately from the formula ( |5.12| ) and from the above 
proof of the theorem 1. Namely, from ( p. 12 ) we get: 



<f> 2 (a) ± l^(o) + - R( ai ) 



R(a r . 



dx = (2m =F l)m 



(5.15) 



QED. 



Some remarks are in order. 

First if Fi^((j)) diverged to infinity faster than 4> 2 then every integral of f(F^ 2 ) in (fUl|) 
would contain E'-dependent infinite series not reducing of course to simple values ±Z7T i.e. 
the relation ( 5.14 ) as well as ( |5.15 ) could not be valid any longer. 



Second one can easily check that if Fi^ifi) do not diverge to infinity faster than eft 2 and 
the potentials V±(x,X) are holomorphic then they have to satisfy the assumptions 1° — 5° 
above. 



Therefore we can draw a conclusion that our assumptions about the potentials V± and 
the functions F\ t 2 fit in some way in with the property of V± being shape invariant. But 
as we have checked they are not determined in some unique way by the shape invariance 
condition (|5.§|), (For example, others than the second power rate of growth of Fi^(4>) with 
4> are allowed by (p\g|), see Appendix 3). 



6 Discussion and conclusions 



In this thesis we have demonstrated that there are two basic symmetries, a reflection: 
x — > —x and a translation: x — » x + a, of potentials and of their corresponding Stokes graphs 
which decide whether the JWKB quantization formulae are only approximations to the exact 
formulae ( |2.iq ) or they are exact by themselves. 

We have established also that despite the above two symmetries for the latter case to 
happen an additional property of the considered potentials and the corresponding Stokes 
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graphs has to be present. Namely, this is the simplicity of SG's generated by the original 
potential expressing itself in no more than two turning points and in no more than two singular 
points in the basic period strips to appear. In the opposite case a proliferation of additional 
sectors in the basic period strip prevents the periodicity properties of the corresponding 
quantization conditions ( 2.1 0| ) to be used to reduce the conditions to the pure JWKB ones. 



The possible relaxation of these conditions has been described in Sec.||. and the corresponding 
examples were given in [|J. 

Altogether, the above two symmetries and the simplicity condition reduce effectively a 
number of exactly JWKB-quantized potentials to only eight of them. All of them have long 
been known. But due to our investigations we have given them the property of being rather 
exceptional. 

We have also shown that the SUSY JWKB exact quantization formulae seem to be only 
different formulations of the exact conventional ones at least in the case of the translationally 
shape invariant potentials. We have supported the validity of this conclusion showing the 
exactness of the JWKB formulae for the two cases of the shape invariant potentials (Vg and 
V\o of Sec.[54]) not found by our earlier analysis of Sec|| Additionally, our two theorems of 
Sec. |5.1| suggest also that there is close relation between the property of being translationally 
shape invariant and the meromorphic structure of the considered potentials on the x-plane 
which is constrained to contain a limited number of second order poles (in the whole x-plane 
or in the basic period strip if the potential is periodic) and to have a particular behaviour at 
the infinity. 

We would like also to stress that the earlier proofs of the exactness of some JWKB 
quantization formulae as done by Rosenzweig and Krieger [Q, [5] and in the case of their 
SUSY forms by Crescimanno [111]] are incorrect by erroneous calculations of necessary phases. 



We have to note also that the results obtained by Inomata et al |27| for the form of the 
SUSY JWKB formulae in the cases of the broken supersymmetric potentials do not contradict 
ours since the latter concern their exact, not approximated forms which appear to coincide 
rather with those of Comtet et al [Ml- 
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Appendix 1 



Here we show that for the following holomorphic 27ri-periodic function: 



q (x,E,\) = Y J q n (E,\y 

n=l 



(Al.l) 



with even k-l and having only simple zeros its Weierstrass product representation is the 
following: 



q(x,E,X) = Ce^x ]J 



1-^ 

Xn. 



(A1.2) 



where C = q(x, E, X)/x\ x= o or C = g(0, E,X) if x = is not a root of q(x, E, A). 

The above formula follows from the observation that Q(x, E, X) = q(x, E, X) ■ 
■ exp(—k/2 — 1/2) is also 27ri-periodic and holomorphic with the same roots as q(x, E, X) and 
therefore its WP representation should be: 



Q(x,E,X) = Ce ax x J] 

n>l 



X 
J i, 



(A1.3) 



where a is an integer by periodicity of Q. On the other hand the representation ( Al.3| ) 



depends analytically on the coefficients q n of (Al.l) and we can always choose them in such a 
way to make Q symmetric under the reflection: x — > —x. This operation does not change in 
( A1.3| ) the product itself (the distribution of roots are then invariant under the operation) but 
changes e ax into e~ ax . However, a being integer cannot change with analytic continuation 
od g n 's and therefore it has to be zero from the very beginning. 



As an example consider q(x, E, A) given by fl3.4|) for which its distribution of roots is 
shown in Fig. 3. We have for it: 



ae 2x - 2pe x + 7 = (a - 2(5 + i)e x J] 



n>l 



X 
Xn 



(A1.4) 



We want to calculate with the help of ( |A1.4 ) a change of phase of q(x, E, X) when trans- 
porting it from a point xq of the line Qx = ir to the point xq — 2iri of the line Qx = — ir. To 
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this goal we note that as it follows from ( A1.4j ) the roots of q(x, E, A) lying in large distances 



from the points considered almost do not contribute to the values of q{x, E, A) in the consid- 
ered strip (their product in (|A1.4 ) is close to 1). Therefore we can take a sufficiently large 



but finite number of roots around the considered points to perform the calculations needed 
(eventually we can take the limit of the infinite number of roots). 

Starting from the point xq we can consider n pairs of roots lying above the line = ir 
(n is large) and n pairs of roots lying below the line. The arguments of xq — Xk we take to 
be positive for Xk lying below the line Qx = tt and negative in the opposite case. It is clear 



that the net result of summing the corresponding arguments of the product in ( A1.4 ) is zero. 
But there is still non zero contribution to the argument of q(xQ, E, A) coming from the factor 
e x of QA1.4 ). It amounts of course to ir and this is the total argument of q(xo,E, A). 



At the points xq — vri our calculations are similar. Keeping the same set of roots as 
chosen previously we see that to the total phase of the product at xq — iri contribute only the 
two most distant pairs of roots lying above the line 9x = tt so according to our convention 
this contribution amounts to 4(— vr/2) = —2tt (in the limit of the root number going to 
infinity). Together with the argument — tt provided by the factor e x we get the argument of 
q(xQ — ici, E, A) to be equal to — 3ir. Therefore the total change of the argument of q(x, E, A) 
between the lines considered is equal to — 4ir. 



Appendix 2 



We demonstrate here particularities of our statement that the exactness of the conven- 
tional JWKB formula for the potential Vg(x, A) = Vg-(x, A), the first one of those in (|5.9|), 
follows from its SJWKB one. For the potential Viq(x, A) our considerations would be similar. 
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Fig. 16 The two sheeted <^>g-Riemann sur- 
face for the unbroken superpotential 4>g 
(the case 1°) 



To this end let us consider the relation ( 5.12| ) using the superpotentials <pg(x,X) given 
above. First consider the case 1° of the exact supersymmetry. The corresponding Riemann 
surface Rj, g is depicted on Fig. 16. This is two sheeted surface with the branch points at 
4> 9 = ±i(a 2 -b 2 )2. The latter are the unique singularities of the integrand of the following 
integral: 



K 



-F x ( 



E 



E 



F 2 (fo 



(A2.1) 



(since roots of F 2 at <f> = ±2£J are also the roots of F\). 



Rfa is, as it can be easily noticed, a map of the basic period strip — it < x < tt of the x- 
plane (see Fig. 15), so that the four turning points of qg-(x, A, E) from this strip are mapped 
pairwise into R^: the two from the segment (— 7r/2,7r/2) into the sheet a) of Fig. 16 and the 
other two into the second one. It is also easy to note that in the quantization formulae Q5,l| ) 
and (|5.2| ) the contour K\ on Fig. 15 can be substituted by the contour K 2 of the figure by 
the periodicity. The contours are mapped into R^ as and respectively, the latter 
surrounding the respective pairs of the turning point pictures on R^ (see Fig. 16). Therefore 
for the quantization formulae (|5.1[) and (|5.2|) we can write: 




■^09 + -Tp ~ Edx 



A 
2 



+ 

K x K 2 



S_ 

A 1 
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Ki 



Edx - - 



A'i 



+ 



A' 



2,0 



TV 



E 



F 2 (<fo) 



Ki 



Edx - - 



K. 



E 



E 



F 2 (fa 



where K oox ^ and K OQ2 ^ are the contours obtained by an obvious deformations of the contours 
K\ 6 and K 2 ^ which contain all the singularities of ^1,2(^9)- Making use of the explicite forms 
of ^1,2(^9) as given in Sec.|5] we can calculate the last integral in ( A2.2| ) getting for it the 
value Altogether with (|5.2|) this gives the result (|5.1|). 



Consider now the broken case 2° of the superpotential 4>g. The corresponding basic period 
strip of qg-(x, X, E) and the quantization contours K\ and K 2 transform into as it is 
shown in Fig. 17. Once again we can write the sequence analogous to ( A2.2j ) deforming the 



contorus ki,<p and ^2,0 of Fig.17 into and respectively to perform the final 

integration getting again +iir and consequently the exact formula (5.1). Of course, the 
starting value of m can be now zero. 



A lm 





Fig.17 The two sheeted 0g-Riemann surface for broken superpotential cpg (the case 2°) 



Appendix 3 

We shall show here that the shape invariance condition (^^) does not prevent in some 
obvious way for F\^.{4>) to diverge with any power of <\> when (ft — ► 00. To this end let us 
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rewrite Q5.S| ) in terms of superpotentials. We get: 



i> 2 (x, A, a) + —<j/(x, A, a) = <fi 2 (x, A, a\) — r<t>'{x, A, a±) + R(a\) 



(A3.1) 



ai = /(a) 



Introducing farther to (|A3.1 ) the function F2(<fi,a) (= (f>'(x(<j),a),a)) we obtain: 



2<f>A(<p,a) + A 2 (<P,a) + R(f(a)) 
2 + A' (6, a) 



(A3.2) 



where A(d>, a) is defined as: 



A(4>, a) = A(x((j), a), a) 



(A3.3) 



(x, ai) = 0(ic, a) + A(x, a) 



It follows from ( A3,2j ) that the behaviour of F^^, a) when <p — ► cxd comes out from the 
corresponding behaviour of A(0, a). The latter, however, under the assumption that <j> = oo 
is at most a pole for it has to be following: 



fc>0 



-fc+i 



(A3.4) 



i.e. this pole has to be simple at most. 



The last equation is a conclusion of the condition: 



x((f) + A(4>, a), ai) = x((f>, a) 



(A3.5) 



under which the shape invariance property ( A3.1 ) is satisfied. Note also that due to equality: 
x'J(j),a) = 1/F2((f>,a), the following relation is coming out from ( [A3. 5 ): 



F 2 {x{4> + A(<j>, a), f(a)) = (1 + A'^, a))F 2 (<t>, a) 



(A3.6) 
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It is now easy to conclude from ( |A3.2 ) that if 6q 7^ 0, —2 then F2(4>,a) grows as (j) 2 when 



— ► oo. But for example if bo = —2 and b%, b n ^0 with 62, ■ •■> &n-i = 0, n > 2, then ^(0, a) 
has to grow as (f> n+l when <f> —* 00. 

Of course, whether A(0,a) can really behave in the above ways depends totally on the 
properties of the superpotentials considered which on their own are constrained by flA3,6j ). 
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